On Retardation Effects in Space Charge Calculations Of High Current Electron Beams 
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Laser-plasma accelerators are expected to deliver electron bunches with high space charge fields. 
Several recent publications have addressed the impact of space charge effects on such bunches after 
the extraction into vacuum. Artifacts due to the approximation of retardation effects are addressed, 
which are typically either neglected or approximated. We discuss a much more appropriate cal- 
culation for the case of laser wakefield acceleration with negligible retardation artifacts due to the 
calculation performed in the mean rest frame. This presented calculation approach also aims at a 
validation of other simulation approaches. 



I. INTRODUCTION 

Laser acceleration is a promising field in various as- 
pects, such as the miniaturization of the accelerator setup 
and the availability of electron beams with high-current 
density. The plasma wake field acceleration in the bubble 
regime was predicted in PIC simulations POl and led to 
rapid progress in various experiments 0, H EL HI • 

The effect of the charge density on the electron beam 
dynamics after the extraction into vacuum can be esti- 
mated when considering the electromagnetic field energy 
@ per particle in the mean rest frame of the particle 
bunch 



1 



4ne N 



N i-1 

EE- 

i=2 j=l I 



(1) 

where N is the number of charged particles considered 
and qi is the charge of the particle i at the position r',. 
Within the mean rest frame, the average electron veloc- 
ity in the longitudinal direction is at a minimum. For a 
large N, a homogeneously charged sphere with a radius 
R scales as u'* oc Q 2 /R with Q being the total charge. 
Laser accelerated electron beams yield normalized field 
energies u n = u'j /m e c 2 which enters a regime which is 
far above the one that can be reached by conventional ac- 
celerators. Space charge effects in this regime have been 
examined under various aspects, these are, for example, 
longitudinal wakefields @, @] , energy spreads introduced 
at the electron extraction from the plasma 0] and the 
temporal development of the induced energy chirp fl(ij |. 

A frequently used method to calculate space charge 
effects consists of simulations based on point-to-point in- 
teractions (PPI). Numerical calculations based on PPI 
are particularly exposed to artifacts if the calculation is 
performed in the laboratory frame. This issue is dis- 
cussed in [ll| explaining the artifacts in the way PPI 
simulations commonly account for retardation: In the 
absence of the knowledge on the particles 4D trajectory, 
retardation is approximated by assuming constant ve- 



* Corresponding author; Electronic 

stefan.becker@physik.uni-muenchen.de 



address: 



locities. Consequently, the acceleration of the simulated 
particles, i.e. the macro particles during each finite time 
step, is neglected. For the regime discussed here, these 
artifacts are examined in detail. Much more appropriate 
results can be obtained when the calculation is carried 
out in the mean rest frame of the electron bunch using 
PPI. In this case, the average velocities of the macro par- 
ticles are only weakly relativistic and consequently the 
retardation artifacts are minimized. This result is then 
discussed in comparison with other widely applied cal- 
culation approaches, which exhibit considerable artifacts 
which manifest in different characteristics of the longitu- 
dinal phase space. Therefore, the design of experiments 
using laser accelerated electron beams could possibly be 
misguided, especially table-top Free-Electron Lasers [HJ 
which crucially depend on the characteristics of the longi- 
tudinal phase space. The calculations here are performed 
using GPT [l3j. The artifacts addressed, however, are 
not specific to the GPT code, but can be found in any 
code utilizing the PPI model. 



II. COULOMB EXPANSION 

The electron bunch considered here is cold and thus has 
zero divergence and zero emittance. The initial bunch 
configuration with a spatial Gaussian density distribution 
in the laboratory frame has the RMS values a x — o y — 
a z = 1/im, an initial kinetic energy given by the Lorentz 
factor 70 = 300 and a total charge of Q — InC and hence 
results in a normalized field energy of u n » 10%. 

For obtaining a 4D trajectory, we have to make as- 
sumptions for the acceleration process: Within the 
plasma, any expansion driven by the Coulomb forces is 
suppressed due to the strong plasma fields. We can hence 
assume that the bunch is in a Gaussian shape before leav- 
ing the plasma accelerator in both relevant frames of ref- 
erence, laboratory frame and the mean rest frame. The 
end of the plasma is assumed to be a "sharp" edge be- 
hind which the expanding effect of the Coulomb forces 
is suddenly switched on. This means that only electrons 
that have crossed the boundary take effectively part in 
the Coulomb interaction. 
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A. Mean Rest Frame 

The mean rest frame for the case considered here is an 
inertial frame of reference co-propagating with the elec- 
tron bunch at the constant normalized velocity (3 a . Be- 
fore the Coulomb interaction is switched on, the Lorentz 
transformation leads to a bunch prolongation of a factor 
of 7o- In this frame of reference, the Coulomb inter- 
action cannot set in instantaneously and globally. In- 
stead, the onset of the interaction starts at the front end 
of the bunch and spreads towards the rear end. Note 
that the bunch geometry in the mean rest frame leads to 
a transverse expansion dominating over the longitudinal 
debunching. 




FIG. 1: The electron bunch is shown schematically in (a) at 
the end of the wakefield acceleration in its initial configuration 
using a Minkowski diagram in the laboratory frame (red) and 
the mean rest frame (blue). The spatial axis z is the longitu- 
dinal propagation direction. The plasma boundary is located 
in the laboratory frame at z = and in the mean rest frame 
at z' — 0. pi, p2 and P3 are 4D trajectories corresponding to 
the rear end, the center and the front end of the bunch, (b) 
shows the beam state using a PPI simulation performed in the 
mean rest frame in correspondence with (a) . The longitudinal 
bunch position s' is plotted against the normalized transverse 
velocity j3' x with s' = being the center of mass, (c) shows 
the situation when half of the particles have left the plasma 
and (d) correspondingly in the mean rest frame at t 1 = 0. 

Fig. [T] shows Minkowski diagrams and particle distri- 
butions in the mean rest frame of the bunch. Panel [TJi 
illustrates the initial beam condition in the laboratory 
frame and in the mean rest frame of the electron bunch 
before leaving the wakeheld accelerator. The bunch is as- 
sumed to exhibit a Gaussian spatial density distribution 
in both frames of reference. This state in the mean rest 
frame is also shown in panel |T|b using a PPI calculation. 

Panel [TJ; illustrates the beam in the laboratory frame 
and in the mean rest frame of the electron bunch at the 
time t = t' = 0, where the plasma boundary is at the 
center of the bunch and moves at the constant velocity 



(3 p b = — (3q. This center of the bunch is also determined 
to be at the longitudinal position zero in both frames 
z = z' = 0. 

In general we can state that the calculation performed 
in the laboratory frame does not require sophisticated as- 
sumptions concerning the initial simulation conditions, 
since the predominant transverse space charge interac- 
tions can be assumed to set in instantaneously: The du- 
ration of propagation of the plasma boundary in the lab- 
oratory frame is ultra fast, which means that the shape of 
the bunch does not significantly change during the prop- 
agation of the boundary. Therefore, simulations yield 
same results for the cases of the electrostatic interac- 
tion being switched on instantaneously and the moving of 
the boundary through the bunch. Within the mean rest 
frame, however, the different ways of coincidence causes a 
longitudinal bunch prolongation, which can be seen from 
panels QJ,,c. Panel QJl shows the spatial electron density 
distribution being altered while the plasma boundary is 
still within the bunch. As a precondition, however, we 
can assume that the bunch has the same spatial sym- 
metry at the end of the wakefield acceleration in both 
frames of reference, i.e. the mean rest frame and the lab- 
oratory frame as shown in panelQJu Since we have chosen 
a Gaussian density distribution, we have a spatial point 
symmetry for the case discussed here. The end of the ac- 
celeration distance is at z — 0, which requires the intro- 
duction of a constraint whereby the onset of the Coulomb 
interaction propagates from the front end to the rear end 
of the bunch as shown in panel [T];. 

The further evolution of the electron bunch calculated 
in the mean rest frame is shown in the laboratory frame 
in Fig. [2 The Lorentz transformation into the labora- 
tory frame is performed assuming constant velocity and 
is given by 

P xl = 7o(l+A)/3V) Xi = X 'i 

= 7o(l+/3o^) Vi = ( 2 ) 
a §U _ A 

Pzi- l+Pop> xi Z i~ 7o(l+A)/3V)' 

where Zi{z'^) is obtained using linear extrapolation in the 
(z, ct) space. A slice is a longitudinal subsection of the 
bunch. The small slice energy spread can be explained 
regarding the longitudinal phase space in the mean rest 
frame (paneldJa) : Particles with highest values of 7' origi- 
nate from the regions with the highest initial electrostatic 
fields, i.e. in the longitudinal direction along the axis 
around the center of the bunch and transversely from off- 
axis. In contrast to the lab frame, the bunch shape in the 
mean rest frame changes significantly already during the 
propagation of the plasma boundary (panel [TJi). The off- 
axis electrons close to the boundary are predominantly 
accelerated towards the rear (left) end of the bunch and 
therefore pushing on-axis electrons towards the head as 
a result of momentum conservation. Due to the Lorentz 
transformation 7 = 7o7'(l+/^) for a particle at a certain 
position s', a larger value of 7' is reduced by a negative 
value of j3' z and a smaller 7' is boosted by a positive (3' z . 
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FIG. 2: The further evolution of the bunch in Fig.[T]is shown 
3 ps later in the longitudinal phase space with the calcula- 
tions performed in the mean rest frame, (a) displays the the 
longitudinal phase space of the bunch transformed to the lab- 
oratory frame using Eqs. The slice energy spread yields 
significantly smaller values than the total energy spread of 
the bunch. This fact also leads to a gap in the phase space. 
The reason for this property can be seen in (b), where the 
longitudinal phase space is drawn with the longitudinal nor- 
malized velocity color coded in the mean rest frame. The 
plasma boundary leads to an asymmetric particle distribu- 
tion in the mean rest frame. The instantaneous onset of the 
Coulomb interaction in the mean rest frame would lead to a 
symmetric particle distribution as shown in (c) and (d) in a 
direct comparison after 0.7 ps of expansion. 



Finally, particles with different values of 7' obtain virtu- 
ally the same 7 in the laboratory frame, which leads to 
the "gap", i.e. the small slice energy spread for electrons 
at a specific bunch position s as can be seen in panel [5^,. 

A symmetric particle distribution in the mean rest 
frame would be obtained from the instantaneous onset 
of the Coulomb interaction and is shown in panels |Sfc,d: 
The phase space does not yield the correlation described 
above and thus, the slice energy spread is larger. 



III. COMPARISON OF CALCULATION 
APPROACHES 

The calculation of the space charge driven expansion 
is examined using PPI according to [y, [l3j, LL4[ . The elec- 
tromagnetic fields are calculated relativistically, where 
radiation effects are neglected and retardation is treated 
in accordance with the constant velocity approximation. 
The Coulomb field of particle j acting on i is given in the 
rest frame of j by 



Qr'i 
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(3) 



i + -±r (rji • 03 % =r'i~ r'j , (4) 



with Q being the charge of the macro particles and r'ji 
being the distance between the particles in the rest frame 
of j. The Lorentz transformation of the electromagnetic 
fields of particle j acting on i in the laboratory frame 
yields 



IjPj x E 'j 



7J -> -> -. 



7i + l 



(5) 



A tracking code [13| applying Eqs. ([5]) is used to cal- 
culate the free drift of the considered electron bunch in 
vacuum. Fig. [3] compares results of different calculation 
methods at a later point in time than Fig. [2] The ap- 
propriate calculation using PPI performed in the mean 
rest frame is shown in the column [3j\ The results being 
Lorentz transformed are shown in column [3j. The slice 
energy spread (panel EJc) is larger compared to the one 
shown in panel f^fi because debunching effects on larger 
propagation distances cause a longitudinal phase space 
as displayed in panel \23p developing towards the one as 
in panel Column [3JI shows the PPI calculation per- 
formed in the laboratory frame. The difference in com- 
parison with columnfSJ can be explained with retardation 
artifacts due to the constant velocity approximation as 
described in [TTf| . A further method besides PPI treating 
space charge is the application of Poisson solvers, which 
evaluate the electrostatic space charge field in the mean 
rest frame of the bunch, where magnetic fields are ne- 
glected occurring due to relative velocities. The Lorentz 
transformation into the laboratory frame introduces the 
magnetic fields. In this respect, the Poisson solvers re- 
ferred to solve the equations of motion in the laboratory 
frame. The calculation time using this method linearly 
scales with the number of macro particles and thus al- 
lows many more macro particles to be considered. We 
evaluated a calculation using PPI performed in the mean 
rest frame of the bunch as it would be obtained from a 
Poisson solver by using the fields 



E' 
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B' 



= 0. 



(6) 



We obtained virtually identical results for the exam- 
ple bunch considered here comparing the cases utilizing 
Eqs. (O and utilizing Eqs. © in the mean rest frame. 
Thus, the case examined here yields an appropriate treat- 
ment with respect to the retardation artifacts when uti- 
lizing a Poisson solver. However, the electromagnetic 
field being obtained in the mean rest frame also requires 
the correct consideration of the 4D trajectories of the 
bunch particles as described above. The result of Pois- 
son solvers assuming instantaneous onset of the Coulomb 
interaction leads to a longitudinal phase space (panel[5J;), 
where the slice energy spread within the laboratory frame 
is shown to be overestimated compared to the case of cor- 



4 



II III 




s' [mm] s [urn] s [urn] s [urn] 

FIG. 3: Calculation of the vacuum expansion of the considered electron bunch is shown with the initial state as illustrated 
in Fig. [T] after the propagation distance of 1.8 m corresponding to 6 ns in the laboratory frame or after 20 ps in the mean 
rest frame. The center of mass is located at s = = ^i(t)si(t)/ Y2i 7i(£)- Panels in column I' show the calculation with 
negligible retardation artifacts in the mean rest frame and panels in column I show the Lorentz-transformed results (Eqs. [2} 
in the laboratory frame. Column II shows the PPI calculation performed in the laboratory frame and column III shows the 
calculation result corresponding to a Poisson solver, both with instantaneous onset of the Coulomb interaction. The longitudinal 
spatial inner bunch position is plotted in the mean rest frame (s') and in the laboratory frame (s). Panel I'a shows the line charge 
density distribution, panels I'b,c show the bunch current. Both of these representations are proportional to the longitudinal 
particle density. Row (b) shows the spatial particle density distribution with the transverse coordinate plotted in the mean 
rest frame (a/) and in the laboratory frame (a;). Row (c) shows the longitudinal phase space. The color coding is linear and 
equivalent to the one of panels (b) and (d) in Fig. [T] 



rect initial conditions applied (panel H^). To our knowl- 
edge, the instantaneous onset of the Coulomb interaction 
is assumed among the vast majority of codes which are 
widely applied and which utilize the method of a Poisson 
solver. In addition, more realistic simulation scenarios 
might involve particle bunches having a notable energy 
spread or divergence. In these cases, velocities in the 
mean rest frame might not allow to neglect the magnetic 
fields, where PPI codes using Eqs. © being performed in 
the mean rest frame of the bunch yield the least artifacts. 
Panels Hfc,d and column [3JII show results as obtained 
from a Poisson solver at a later point in time. Coulomb 
interaction is considered in the mean rest frame and is 
assumed to set in instantaneously. This wrong initial 
condition leads to the difference compared to column [3J. 
Considering energy conservation for the example beam 
in Fig. d] yields the total energy being conserved in the 
case of the calculation performed in the mean rest frame 
of the electron bunch and for the case using the Poisson 
solver. 
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FIG. 4: The solid line shows the kinetic energy. The dashed 
line shows the total energy which is the sum of kinetic energy 
and field energy using Eq. ([1]). The calculations correspond 
to Fig. H 



IV. CONCLUSION 

The method discussed allows the calculation of the 
electron bunch evolution in a space charge regime which 
could be reached by laser acceleration with negligible re- 
tardation artifacts. This calculation is performed in the 
mean rest frame of the bunch, where the relative ve- 
locities are only weakly relativistic. Moreover, the re- 
quired Lorentz transformation between the laboratory 
frame and the mean rest frame is shown to be non-trivial, 



5 



since assumptions concerning the 4D trajectories of the 
particles of the bunch have to be made. The result of this 
calculation is compared with the results of two commonly 
applied methods, one using PPI performed in the labo- 
ratory frame and one using a Poisson solver. For both 
approaches, we found significant deviations concerning 
the characteristics of the longitudinal phase space, which 
could ultimately mislead the design of applications using 
laser accelerated electron bunches: A PPI simulation per- 
formed in the laboratory frame principally suffers from 
retardation artifacts, which leads to the violation of en- 
ergy and momentum conservation and to a wrong spatial 
density distribution. The temporal development of the 
energy chirp within the longitudinal phase space is over- 
estimated and thus, the examinations in describe an 
upper boundary. The method of using a Poisson solver 
in principle suffers from a Lorentz transformation which 
has to be done in every time step. Some implementations 
were found to apply the Lorentz transformation for the 
longitudinal spatial position by merely linearly stretch- 
ing the bunch by the Lorentz factor corresponding to 
the velocity of the mean rest frame. Using z(z') as in 



Eq. [2j instead, helps improving the result. This transfor- 
mation, however, assumes constant velocity which might 
introduce artifacts due to the Coulomb driven bunch ex- 
pansion. The method of using a Poisson solver, hence, 
could principally deliver better results than PPI in the 
laboratory frame. However, many simulation codes do 
not allow to include the correct initial conditions. The 
method presented here neither requires Lorentz transfor- 
mations after each time step nor suffers from retarda- 
tion artifacts. Hence, it offers an appropriate calculation 
method for the regime of significant space charge effects 
as expected for laser electron acceleration. 



Acknowledgments 

We thank M. Dohlus (DESY) for fruitful discussions. 
This work has been funded by the DFG through tran- 
sregio TR18 and supported by the DFG Cluster-of- 
Excellence Munich Center for Advanced Photonics MAP. 



[1] A. Pukhov and J. Meyer-ter-Vehn, Appl. Phys. B 74, 355 
(2002). 

[2] S.P.D. Mangles et al, Nature 431, 535 (2004). 

[3] C.G.R. Geddes et al, Nature 431, 538 (2004). 

[4] W.P. Leemans et al, B. Nagler, A.J. Gonsalves, Cs. Toth, 
K. Nakamura, C.G.R. Geddes, E. Esarey, C.B. Schroeder, 
S.M. Hooker, Nature Phys. 2, 696-9 (2006). 

[5] J. Osterhoff et al, Phys. Rev. Lett. Phys. Rev. Lett. 101, 
085002 (2008). 

[6] J.D. Jackson, Classical Electrodynamics John Wiley & 
Sons (1999). 

[7] R.A. Bosch, Phys. Rev. ST Accelerators and Beams 10, 
050701 (2007). 

[8] G. Geloni, E. Saldin, E. Schneidmiller and M. Yurkov, 
Nucl. Instrum. Methods Phys. Res., Sect. A 578, 34 



(2007). 

[9] G. Stupakov and Z. Huang, Phys. Rev. ST Accelerators 
and Beams 11, 014401 (2008). 
[10] F. Griiner, C. Schroeder, AR. Maier, S. Becker and Y.M. 
Mikhailova, Phys. Rev. ST Accel. Beams 12, 020707 
(2009). 

[11] G. Fubiani, J. Qiang, E. Esarey, W.P. Leemans, G. 
Dugan, Phys. Rev. ST Accelerators and Beams 9, 064402 
(2006). 

[12] F. Griiner et al, Appl. Phys. B 86, 431 (2007). 
[13] General Particle Tracer, www.pulsar.nl. 
[14] B. van der Geer and M. de Loos, PhD. thesis, Technische 
Universiteit Eindhoven (2001). 



